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Abstract 

It is shown that for every problem within dimensional regularization, 
using the Integration-By-Parts method, one is able to construct a set 
of master integrals such that each corresponding coefficient function 
is finite in the limit of dimension equal to four. We argue that the 
use of such a basis simplifies and stabilizes the numerical evaluation 
of the master integrals. As an example we explicitly construct the 
e-finite basis for the set of all QED-like four-loop massive tadpoles. 
Using a semi-numerical approach based on Pade approximations we 
evaluate analytically the divergent and numerically the finite part of 
this set of master integrals. 

The calculations confirm the recent results of Schroder and Vuori- 
nen. All the contributions found there by fitting the high precision 
numerical results have been confirmed by direct analytical calculation 
without using any numerical input. 



^On leave from Institute for Nuclear Research of the Russian Academy of Sciences, Moscow, 
117312, Russia. 

^On leave from BLTP, Joint Institute for Nuclear Research, Dubna, Russia. 



1 Introduction 



The evaluation of higher order perturbative corrections plays an important role 
in testing the Standard Model on a high precision level as well as in the search 
for new physics beyond the established theory. Especially given the increasing 
resolution of modern experiments, it is mandatory to push ahead the efforts on 
the theory side. Many reviews, see e.g. [1,2], deal with the application of such 
calculations. The present work focuses on the actual evaluation of the integrals 
appearing when considering higher order corrections. 

Let us begin by noting that basically every analytical method applicable for 
higher order calculations can effectively be put into action only within the "cal- 
culation friendly" environment provided by dimensional regularization [3-5]. Of 
course, computing e.g. the 1-loop f3 function one could use a variety of ap- 
proaches. However already at 2-loop level it is almost impossible to avoid the use 
of dimensional regularization unless one is considering a simple scalar toy model. 

It should also be stressed that from a purely calculational point of view the 
most valuable feature of dimensionally regularized integrals is their similarity 
with the usual absolutely convergent integrals. As a result almost all standard 
operations are mathematically justified. 

In particular, the traditional Integration By Parts (IBP) method is based on 
the possibility to discard the surface terms within this framework^ and allows for 
generating exact identities between different integrals. In a next step the resulting 
equations are used to reduce every Feynman diagram to a limited set of mutually 
independent so-called master integrals. 

This reduction procedure is a highly non-trivial problem by itself as the grow- 
ing number and complexity of IBP identities defies any attempt of a straightfor- 
ward solution starting from some number of loops and external legs. Therefore the 
classical way of "manually" examining the IBP identities seems to be unsuitable 
for such problems. A prominent exception is the crack-down of the long-standing 
problem of computing the generic moments of the structure functions in deep- 
inelastic scattering at three-loop level [9] which has been completed only recently. 
At the moment there are a few approaches trying to automatize the procedure 
of reduction to master integrals [10-22]. Once the reduction has been achieved 
one is left with the evaluation of the resulting master integrals. A complete up- 
to-date review about various methods of calculating Feynman diagrams can be 
found in [23]. 

■^The possibility is by no means obvious for a general case of integrals containing both UV 
and IR divergences and requires a dedicated proof. The latter, in turn, is not possible without 
a proper definition of the very dimensional regularization for such integrals. To the best of our 
knowledge the elaborated mathematically meaningful definition of dimensional regularization 
and its main properties was presented in [6, 7] and reviewed in [8]. 
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As is well-known within the method of dimensional regularization, the usual 
UV divergences manifest themselves as poles in e = (4 — d)/2, with d being 
the spacetime dimension. The calculation of the master integrals is complicated 
by appearance of unphysical, so-called spurious, poles in the coefficient functions 
which are often induced during the reduction procedure. Through the appearance 
of these spurious poles one is forced to compute extra terms in the e-expansion 
of the master integrals. In general the evaluation of every next term in the e- 
expansion is challenging, especially if the application of numerical methods is 
necessary (see, e.g. [24-30]). 

On the other hand, the choice of master integrals is by no means unique. 

In the present work we discuss a convenient way of finding a new basis of 
master integrals. In particular we show in the next two sections that for every 
problem reducible with the help of the IBP method to a set of master integrals the 
latter can be always chosen in such a form that all the corresponding coefficient 
functions are finite in the limit e — > 0. We call this basis the e-finite basis. In 
section 4 we explicitly construct such an e-finite basis for the set of all QED- 
like four-loop massive tadpoles appearing in the first moment of the low energy 
expansion of the vacuum polarization in QCD. Using a semi-numerical approach 
based on Pade approximations we evaluate numerically all the masters of this 
e-finite set and confirm the recent results of Schroder and Vuorinen [26]. Then 
we focus on analytical results which can be easily derived following this method. 
Our general comments are presented in section 5. 

2 Generalities and Notation 

Let us consider a quantity P which we want to compute up to L-loop level. The 
sum of all bare contributions Pg is given by 

PB=i:h'pi'- (1) 

1=0 

Here we have introduced a formal parameter h counting the number of loops. 
The renormalized L— loop contribution is expressible through 

/i^ P(^) = /i^ P(^) + J2 h'Z^'^ h'^-'P^t'^ (2) 
1=1 

as a linear combination of the bare quantities P^^ and the corresponding rcnor- 
malization constants, here summarized in Z^^^ for the /-loop contribution. We 
assume that the calculation of the integrals up to h^'^ order is well under con- 
trol. This includes the sum in the second term of equation (2). Therefore we 
concentrate on the bare L-loop contribution P^\ 
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After the evaluation of all the traces in d dimensions each Z-loop contribution 
can be written as a linear combination of d-dimensional integrals FP(e) and 
associated coefficient functions ff\e) 

i 

Here possible divergencies in e can only result from the evaluation of the various 
integrals Ff\e) whereas the functions ff\€) remain finite in limit e — > 0. The 
resulting number of integrals {F^\e)} is in general extremely high and therefore 
reduced by use of IBP relations to a small set of N master integrals {Mf\e)} of 
Z-loop order. In this framework the contribution p'^^ can be written as 

N 

Pi' -Y.^\e)Mf\e). (4) 

i=l 

Through the dependence on the spacetime dimension d of the IBP relations the 
coefficient functions al\e) in equation (4) can contain divergencies for e — > 0. In 
this case the expansion in e of the corresponding master integrals Mf' (e) has to 
be evaluated deeper than to the finite contribution in order to compensate for 
spurious poles. In the next section we will construct a different set of master 
integrals {Mf\t)} such that the expansion 

TV 

P« = $:c«(e)Mf(e) (5) 

i=l 

leads to coefficient functions c'Pie) being finite in the limit e ^ 0. This new set 
of master integrals is called e-finite basis. 



3 Construction of the e-finite set 

The comparison of equations (3) and (4) 

N 

suggests that the e-finite basis can be found in the set of integrals {Ff\e)} as 
the coefficients ff\e) already expose the desired behavior. By definition all the 
initial integrals can be expressed in terms of master integrals 

N 

Ff\^) = Y.^h^)Mf\e). (7) 
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Among this in general huge set of equations only the subset {-^di'v} exhibiting 
a divergent coefficient function a^lj{e) is interesting. Let amn(e) be one of the 
functions showing the l/e"™*^ -divergency with the highest power nmax- Starting 
from the relation for Fm{e) involving amn(e) we are able to construct a new basis 
of master integrals by solving equation (7) for the master integral associated with 



(8) 



As amn(e) was divergent for e — > the reciprocal function is finite. By means 
of this replacement Fm{^) is now a member of the new set of master integrals 
{Mf\e)} replacing M^\e) in the following way: 

{Mf\e)} = {M«(6), . . .,Mt,{e\F!f^{e\Mil,{e\ . . . ,M«(6)}. (9) 

In a next step the replacement (8) is performed for all equations in the set {-^^divl- 
This procedure is iterated until all coefficient functions in {-Fdivl finite. 
The thus obtained final set of master integrals {Mf\e)} fulffils the relation 

Pi' = Y.f?^p^^ = (^)' (10) 

where each appearing coefficient function is by construction finite. If we now 
identify 



(0 



i 

MP{e) = MP{e), (12) 

we arrive at the desired form (5) and call the set of integrals {M-^\e)} the e-finite 
basis. We again stress the fact that all the integrals in this e-finite basis need to 
be evaluated only up to the finite part in their e-expansion, as by construction 
no spurious poles arise. 



4 Applications to four-loop tadpoles 

The concept of the e-finitc basis is now applied to the set of master integrals 
appearing in the calculation of QED-like four-loop tadpole diagrams. We con- 
centrate on the subset of master integrals shown in figures 1 and 2 needed for 
the calculation of the low energy expansion of the four-loop vacuum polarization 
function, completing the work in [34]. We denote each of the master integrals 
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41 



T51 




Figure 1: Subset of master integrals fully analytically known and therefore not 
considered in the construction of the e-finite basis. 



T52 

:-2,4) 




Figure 2: Subset of master integrals for which the e-finite basis is constructed. 
The two numbers in brackets {ni,n2) are decoded as follows: rii is the maximal 
power of a spurious pole (^)"^ which could appear in front of the diagrams pic- 
tured above; n2 is the maximal analytically known power e"^ of the e-expansion 
of the same integral as determined in [31,32] (for T52) and in [26,33] for all the 
rest. 



according to the following rule: after a capital letter "T" we write the number of 
lines in the given diagram. The second number enumerates the different topolo- 
gies exhibiting the same number of lines. The diagrams of the first subset shown 
in figure 1 are analytically fully known and expressible in terms of F-functions. 
We thus construct the e-finite basis only for the second subset of the most compli- 
cated master integrals shown in figure 2. This is done by treating the coefficient 
functions of the "simple" diagrams in the set of equations {-F^fvl finite. 



4.1 Construction of the e-finite set 

Following the description in section 3 we find the diagrams shown in figure 3 as 
the e-finite basis. Here the name of the e-finite master integrals is deduced by 
using the name of the original master, which had been replaced, and the addi- 
tional superscript letter "f " . Note that except for T^^ all members of the e-finite 
basis pictured in figure 3 are just scalar Feynman integrals. The integrand corre- 
sponding to T/4 contains a scalar product of two line momenta in its numerator. 
According to equation (7) the integrals of the e-finite basis can be expressed in 
terms of the original master integrals of figure 1 and 2. This has been achieved 
by using the traditional IBP method in combination with the Laporta Algorithm. 
In this context the IBP identities have been generated with numerical values for 
the powers of the propagators and the irreducible scalar products. In the next 
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Figure 3: Diagrams of the e-finite basis. 

step, the resulting system of hnear equations has been solved, by expressing sys- 
tematically complicated integrals in terms of simpler ones and substituting the 
solutions into all the other equations. For this reduction to the original master 
integrals an automated program, which is based on FORMS [35, 36] has been de- 
veloped in which also partially ideas described in ref. [11,17,37,38] have been 
implemented. The coefficient functions of the integrals in the linear system of 
equations depend on the spacetime dimension d and need to be simplified during 
the solving procedure. This has been done with the program FERMAT [39]. Having 
expressed the integrals of figure 3 in terms of the original basis and following than 
equation (8), one obtains the relations which replace the original master integrals 
by the e-finite ones. As these equations are rather lengthy, they are summarized 
in appendix A. 

4.2 Calculation of the e-finite basis 

In order to calculate the results for the e-finite basis we follow the method based 
on the ideas in ref. [40] and already described in [41]. One line of the four-loop 
tadpole under consideration is cut, thus obtaining a three-loop self-energy dia- 
gram. Since there is no other four-momentum involved, the self-energy diagram 
only depends on the momentum flowing through the applied cut. Therefore we are 
able to immediately perform the angular part of the remaining loop integration. 
Following this procedure, the calculation of the four-loop tadpole is transformed 
to the problem of evaluating a three-loop self-energy and a one-dimensional nu- 
merical integration over the modulus of the momentum through the cut. 

The resulting three-loop diagrams are treated with the help of asymptotic 
expansions. Using the program packages EXP [42,43], MAT AD [44] and MINCER 
[45] one is able to compute the low and high energy expansion of the self-energy 
diagram. In a final step we combine the obtained high and low energy expansions 
in a Fade approximation as described in [46] . The resulting function is integrated 
numerically up to a cutoff A as the integration up to infinity is in general diverging. 
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The remaining integration region from A to infinity can be treated by direct 
analytical integration of the high energy expansion. Therefore one is able to 
compute the structure of the e-poles analytically. We checked that the result for 
the finite part of the master integrals is cutoff- independent. In addition we tested 
that for finite tadpole diagrams like Tg^ the difference between the introduction of 
A and the numerical integration up to infinity is below the error of the numerical 
integration. A more detailed description of the procedure can be found in [47]. 
For diagram Tg^ we explicitly show the steps in the next subsection. 

Note that in general we have a vast freedom in choosing the members of the 
e-finite basis. Thus we give some additional remarks about selecting appropriate 
diagrams. For the applicability of the Fade method in the described way the self- 
energy diagram resulting after cutting one line should not contain any massless 
cuts (see [47]). Nevertheless, for each integral shown in figure 3 at least two 
possible cuts leading to two different self-energy diagrams could easily be found. 
This also provides us with a additional possibility to check the results. 

The self-energy integrals appearing after cutting one line might contain UV 
and/or IR divergences. In general this is not a problem for the presented method 
as one is always able to subtract suitable diagrams containing the same sub- 
divergencies but of lower loop order, thus ending up with a self-energy function 
finite in the limit e — 0. In appendix B we demonstrate this procedure explicitly 
for one diagram. However, a simple power counting applied to the diagrams 
of figure 3 demonstrates that excluding diagram T^^ every possible cut leads to 
finite self-energies. In fact, the minimal presence of any divergencies was one of 
the main principles used by us in choosing suitable candidates for the e-finite 
basis. 

4.3 Treatment of the basis integral T^2 

This diagram has two possible equivalent cuts which lead to the same finite self- 
energy function S"^52 (p) related to the diagram 



resulting from the cut procedure. In general one might be left with a diagram 
which does not allow for a cut resulting in a finite self-energy function. In ap- 
pendix B we present an alternative proceeding for diagram T/2 involving a cut 
leading to a divergent self-energy function. As outlined in the previous section 
we calculate the high and low energy expansion using EXP, MATAD and MINCER 
before combining these results in a Fade approximation 




(13) 



ao + aiz + ■ ■ ■ + ttnZ' 



.n 




biZ -\ h bmZ' 



(14) 
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Given n\ow terms of the low energy and rihigh terms of the high energy expansion 

rpf 

we are able to evaluate a Pade approximant '^[n]m]iP ) self-energy function 

in equation (13) for all values of n and m obeying the relation 

riiow + rihigh = n + m. (15) 

This function we are then able to integrate numerically up to a cutoff A according 
to 

-iWl '"TT^^^^-MiP') (17) 

where mcut = m if we cut a massive and mcut = if we cut a massless line. 

The remaining integration from the cutoff A to infinity is performed by explicit 



analytical integration over the high energy expansion Sp3oo(p^) of the self-energy 
(13). In contrast to (17) here one just has to replace the Pade approximant by 
the high energy expansion 

9 d/2 roo d-1 



For mcut = or if we expand the denominator of the integrand in case of mcut = "t. 
this integral has the simple form 

^max /»00 

4,,^oo(A) = E/ dpa,{e)p--^^^^ (19) 
1=0 -^^ 

with noninteger positive exponents ni{e), at most linearly dependent on e, and 
can be evaluated analytically 



■'111 CIA / \ 

r/2„_(A) = E z;^!^ (A)--^^)^^ (20) 



As some of the exponents ni{e) are linear dependent on e, one has to evaluate 



the high-energy expansion Sp3oo(p^) up to the order linear in e. Following this 
way the divergent part can be fully recovered through this analytical integration 
as we have chosen the diagrams and the cut in such a way that the self-energy is 
finite for e — > 0. The final result for the diagram T/g is then given by the sum 

^5^2[n/m] ^ ^54[n/m](A) + ^54,g-*oo(A) (21) 

of the two contributions (17) and (18). Here numerical evaluations only enter in 
the calculation of the first term resulting from the Pade approximation. Therefore 
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J? 

52niow 




7 


?^high 
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4 


27.30068104(1290) 


27.30068323(70) 


27.30068178(198) 




5 


27.30067759(59) 


27.30067456(1180) 


27.30068017(86) 


6 


27.30068220(998) 


27.30068156(146) 


27.30068059(36) 




7 


27.30067826(492) 


27.30068069(34) 


27.30068067(5) 



Table 1: Results for diagram T52 depending on different number of low and high 
energy input terms. In brackets we give the error of the procedure deduced from 
the standard deviation. 



we obtain the full pole part of the diagrams analytically. In addition, we stress 
again, that the evaluation of any higher term in the e-expansion is not needed. 

According to equation (15) we are left with some freedom how to choose the 
Fade parameters n and m for a given value of niow and rihigh- Thus we average 
for each combination of niow and Uhigh the results of all possible n, m- values. The 
average value tL^ ^ is what we give as the final numerical result for a given 
set of high and low energy terms. The resulting standard deviation is a measure 
of the accuracy of this numerical approach. The values obtained for diagram 
depending on different values of riiow and riMgh are given in table 1. The result is 
nicely converging for increasing number of input terms. 



4.4 Results for the e-finite basis 

Following this method we evaluated all the diagrams of the e-finite basis shown 
in figure 3 up to the finite part. The numerical value of the constant term is 
obtained by using again nine terms in the high and seven terms in the low energy 
expansion of the self energy function arising after the cut. For each integration 
we use for the euchdian integrals in analogy to [34] 

~irjj2^--- ^^^^ 

with the normalization factor 
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In the following we consider m — 1. The results for the e-finite basis then read: 

f _3C3 15C3 91 tt" 47rMog^2 41og^2 
2e2 2e 360e 3e 3e 





32Li4(^) 
+ + 27.30068067(5), 


(24) 


-'54 


^^^^ 6.762240238547(2), 


(25) 


61 


5 

= 29.703462427815(3), 


(26) 


-^62 


= 2.440345823350757(6), 


(27) 


-^64 


= ^ 18.026245978729184(4) , 


(28) 


-^71 


= 5.19831391(6), 


(29) 


-^72 


= 26.4964794044474738(5), 


(30) 


-'si 


= 22.53760348(3), 


(31) 


rrf 
-'-91 


= 1.808879546207(2). 


(32) 



We emphasize that all divergent parts of these integrals are obtained analytically 
without any numerical input and that higher orders in the e-expansion are not 
needed. A further investigation on the analytic results presented here is given 
in the next section. By applying the relations from appendix A the results for 
the e-finite basis can be related to the master integrals in figure 2. We find full 
agreement with the values found in [26] as it is shown in table 2. 

4.5 Analytical results for the master integrals 

In subsection 4.3 we discussed a method which can be used to find the pole part 
of any massive four-loop tadpole analytically. More generally, this procedure is 
formulated in the following statement: the pole part of any (L -|- l)-loop tadpole 
is obtained analytically in terms of pole and finite parts of properly constructed 
L-loop tadpoles and L-loop massless integrals. An explicit construction of such 
an algorithm for the particular case of three-loop completely massive tadpoles 
can be found in [48]. 

It is instructive to investigate in some detail the internal mechanism of the 
cancellation of the spurious poles after the "old'' basis set is expanded in terms 
of the new e-finite basis. As we will see, the very existence of such cancellations 
leads to additional analytical results for master integrals. 

With Sn we denote the maximal power of spurious poles which could appear 
in the coefficient function of an "old" master integral T^. The values of for 
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diagram 


Pade method 


rr-\/^l 

[26] 




27.30068067(5) 


27.30068068 


rpf 

-'54 


-6.762240238547(2) 


-6.762240238546 


mf 


-29.703462427815(3) 


-29.703462427816 


rpf 

-^62 


2.440345823350757(6) 


2.440345823350754 


rpf 


-18.026245978729184(4) 


-18.026245978729184 


Tk 


-5.19831391(6) 


-5.19831389 




-26.1961791011171738(5) 


-26.1961791011171733 


ryf 

-'si 


-22.53760348(3) 


-22.53760347 


ryf 
-'91 


1.808879546207(2) 


1.808879546208 



Table 2: Comparison of the constant part of the e-finite basis with results ex- 
tracted from [26] . Note that in the last column we only keep the accuracy needed 
for this comparison, whereas in [26] higher precision results are available. 



all non-trivial masters is given in figure 2. Let us start with integral T71. As we 
find a value of Sji = 1 one expects that all coeflicient functions should be at least 
proportional to e when expressing T71 in terms of e-finite basis. Here we, of course, 
exclude terms proportional to the analytically known integrals T41, T51, T53 and 
which had by construction been excluded from the e-finite basis. The explicit 
result reads 

Tn = e (4 T/^ + 6 T,{ + 2 T/J + 0{e^) + . . . , (33) 

in full agreement with our expectations. Here we denote by O(e^) terms with 
coefficient functions suppressed by at least as e^ at e — > and . . . represents 
contributions coming from analytically known master integrals. Equation (34) 
also carries important information about the finite part of T71. Indeed, as the 
pole parts of the new master integrals are known analytically, one could easily use 
equation(34) to get again analytically the pole and the finite part of T71 leading 
to the result 

1 3 26 Cs 41 IOC3 vr^ ^039 
^ "6?"2i^-3i^-^-T-^-60^'^+^^-^-^+^^^^+^^'^- ^^^^ 

Repeating the same procedure for the integrals T52, T54, Tgi, T62, and we 



11 



obtain the following results: 

T52 = - T/2 + 0(e=^) Tf + ... (35) 

_ _ 2 _ 29 _ 163 _ 601 _ 635 _ 112 Cs /4549 136 tt^ 
~~?~3e3~6e2~l27~^ 3~ ^ ^ V^6~ ^ 45 

, 32 7rMog^2 32 ^2 952C3\ ^.2^ 
+ o 256Li4(2) ^ +C'(e), (36) 



T54 =e' ^ {t/, + Ti + 1 T/,} + O(e^) + . . . (37) 

1 _ 9 _ 415 4991 _ 89383 /_ 3367679 1792(3 
~ ~ ? ~ 2^ ^ 36e2 ~ 216 e ~ 1296 ^ ^ V 7776 ^ 9 
^2/ 137735095 68992(3 4352 tt^ 10247r2log^2 



46656 27 135 9 

2 8192Li4(i)\ 
+ + :^^^)+0(e^), (38) 



Tei = 1 24 T/2 + 6 T/^ + 12 T/4 + 18 + 6 T/2 } + O (e^) + . . . (39) 
3 19 67 127 3(3 237 _ , vr^ / 3969 



2e4 2e3 2e2 2e e 2 ^ 20 V 2 

559 TT^ 



32 TT^ log^ 2 + 32 log^ 2 + 768 Li4 (|) + 1261 Cs 



60 

+ 189C5)+0(e2), (40) 

= 6^ 1 16 T/2 + 16 T/2 + y T/i I + O(e^) T^' + . . . (41) 

2 4 38 44 16 C3 ,.0 00 > 4^^ / 
3e^ e'^ 3e^ 3e 3e 15 V 



3747r^ 647rMog^2 64 1og^2 2152C; 



+ -^ + 512Li4(|) + 



53 



45 3 

+ 96(5)+ 0{e^), (42) 



T64 = 9e' T/i + 0(e^) + . . . (43) 
1 3 9 2 13C3 147 57r^ / 



27 , , 693C5 



+ 45(3 + ^ +0(6^), (44) 
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= e |2 T/4 + 3 T/, } + 0{e^) + . . . (45) 

_ 1 3 26 Cs 41 Cs 1039 29 vr^ 

~~6e4~2^~372~2^~T^6e~ 120 e 6 120" 

+ 1^^!^^0(.). (46) 

It is remarkable, that depth of the expansion in e of these analytical results for 
"old" masters 754, Tei, ^62, Tq4, T71 and T72 listed above are in full agreement with 
the depth obtained in [26] by fitting the high precision numerical results. 

It is even more remarkable that in mostly all of the cases the highest power e"^" 
of the analytically accessible part of the expansion is correlated with the maximal 
power of the corresponding spurious pole s„ introduced before, such that 

4 = Sn - 1- (47) 

Thus, we observe some rationale as for how many terms of the e-expansion of a 
master integral can be computed analytically. 

However, equation (47) seems to be violated in case of the integrals T54, Tq2 
and T91. Here we observe the relation 

di^Si- 2. (48) 

At first glance equation (48) is in apparent contradiction to the very concept 
of the e-finite basis. Indeed, consider, for example, the integral T54. Given the 
value S54 = 4 from figure 2, it is known that this integral sometimes appears 
accompanied by a spurious pole ^. A naive application of equation (37) tells us 
that after expanding -^T^i in terms of the e-finite basis the integrals T/2, and 
would have non-e-finite coefficients! The problem is solved by noting that 
these ^ poles always appear in one and the same combination, namely 

T = ^ {9 r54 - 48 e + 256 Tgi } . (49) 

As expected the result of expanding T in the e-finite basis has finite coefficients 
in front of non-trivial masters: 

T = 11 r/2 - 8T/4 - 8T/1 + 11 r/2 + 6Ti, + 18r/i + lOr/2 
73 

-^n, + 0{e)Tf + ... (50) 

Hence, this relation should lead to an extra analytical constraint on T52, Tq2 
and T91. Indeed, equations (24-32) allows us to compute analytically the pole 
part of the right hand side of (50) and, simultaneously, the same pole part could 
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expressed through Taylor expansion the right hand side of (49) at e — > 0. As a 
result we get 

It +-^T --T - 5511907345 89 tt^ tt^ 
4 ^^'^ 1024 ^"^'^ 64 ^^'^ 31850496 1440 56 

-ii^ + ii^.44U(.)-i^^-tf-HIii^. (51) 

where we have used the notation: 

oo 

Tn^ (52) 

The integrals T54 and Tq2, being of a "sunset" -like type, are certainly much more 
suitable for a numeric evaluation than T91 (for an exhaustive review see, e.g. 
[49]). Moreover, the integral T91 seems to be rather complicated for a numerical 
treatment as illustrated by the fact that in [26] it was computed with significantly 
less precision than all others. Thus, it is natural to take (51) as an equation for 
the calculation of T91 q. Substituting Tq2,2 and 3 as found in [26] one obtains 
immediately without any further calculation 

Tgi^o = 1.808879546208334741426364595086952090 (53) 

which should be compared to the original result given in [26] (for further discus- 
sion see the next section) 

^91,0 = 1.80887954620833474. (54) 



5 General Comments 

There is a significant amount of freedom in the choice of specific members of the 
e-finite basis. As we have seen in subsection 4.2 this freedom can be used in order 
to greatly simplify the evaluation of the master integrals. In fact, the very idea of 
the e-finite basis appeared after facing problems in trying to evaluate the "old" 
master integrals with the help of the Pade method. 

An important feature of the e-finite basis is extremely simple structure of the 
UV and IR divergences of its members. Indeed, except for Tgi and T91 all members 
of the "old" basis set do suffer from overlapping UV divergences. Whereas all new 
masters but one arc cither finite or possess only mild superficial UV divergences. 
The integral T/2 is not primitively divergent but its only UV divergent subgraph 
could be very easily isolated as demonstrated in appendix B. 

The fact that the IBP identities sometimes provide some nontrivial constraints 
on the values of master integrals is, of course, not entirely new. For example, in 
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[50] the divergent parts of some of the irreducible four-loop completely massive 
tadpoles, appearing in the calculation of the four-loop QCD (3 function, have been 
obtained from the finiteness of the same integral with higher powers of denomi- 
nators. However, we have already discussed in section 4.5 that the evaluation of 
the pole part of a given four-loop tadpole is not a problem at all. We are only 
aware of few cases when IBP identities were actually used to evaluate the finite 
term or terms proportional to a nonnegative integer power of e in the e-expansion 
of nontrivial four-loop master integrals. 

In [33] the order e coefficient of the master integral T54 (there denoted as BB4) 
could be deduced, as it had been the only analytically unknown term entering the 
cancellation of the divergent contributions during the renormalization procedure 
when calculating the so-called p parameter. Furthermore, the coefficient of order 
was subsequently obtained by combining the numerically known value with the 
basis of transcendentals from [51]. Within our approach a// the terms up to the 
contribution are directly accessible by the use of e-finite basis without any need 
for either a high-precision numerical evaluation or the basis of transcendentals 
which could appear. 

Schematically, the relations which were used to fix T54 1 in [33] can be written 

as 

Pb — (known UV sub-divergencies) = 0{e^) (55) 

and 

Pb — (known UV sub-divergencies) = T^i,i/e-\- known terms (56) 

with known UV sub-divergencies. Equations (55) and (56) directly leads to the 
value of T54^i. Nevertheless relations hke (55) involving the cancellation of diver- 
gencies can not, obviously, provide any new information in addition to what is 
obtainable from the use of the e-finite basis. Furthermore there is no systematic 
way to find such relations, not to mention the requirement that all renormaliza- 
tion constants have to be known analytically. 

The fact that such relations can sometimes be found by chance is also well 
illustrated in the next example. In section 4.5 equation (51) was derived in a 
regular way as natural by-product of the properties of the constructed e-finite 
basis. Recently a paper has appeared where the same relation (51) has been 
independently obtained in a similar way like the result for T541 in [33]. The 
authors of reference [52] have computed the so-called decouphng function C,g for 
the QCD quark gluon coupling constant in four-loop approximation. The bare 
function Q is essentially composed of three pieces: 



where the three auxihary bare functions on the right side , C3 ci^id C3 are related 
to the heavy quark containing contributions to the ghost-ghost-gluon vertex, the 
gluon propagator and the ghost propagator respectively (for more details see [33, 
53, 54]). Once the three bare functions are determined and the physically relevant 
quantity is constructed one can renormalize directly with the help of (long- 
available from [55-57]) renormalization constants for and the quark mass and 
without any necessity to renormalize the three auxiliary functions themselves. The 
latter is certainly possible but requires knowledge of three extra renormalization 
constants, namely those for the ghost-ghost-gluon vertex, the gluon and ghost 
wave functions. The corresponding calculation is certainly possible but is not a 
small enterprise and the results for these constants appeared only recently [50, 58]. 
The remarkable fact is that the combination T (defined in eq. (49)) enters into 
the renormalization of the auxiliary functions as follows {b is a known constant) 

Q — (known UV sub-divergencies) does not contain T, (58) 

C3 — (known UV sub-divergencies) = 1 + ba^T + known terms, (59) 

and 

(■3 — (known UV sub-divergencies) = 1 — 2balT + known terms. (60) 
We observe that 

(i) the relation (51) can be found in the process of renormalizing one of eqs. (59 
- 60) but cannot if one deals only with the physical quantity (g. 

(ii) our way of finding (51) is completely self contained and does not require any 
knowledge of extra non-trivial renormalization constants, while the same 
result obtained in [52] does depend on such knowledge. 

(iii) full agreement between relation (53) obtained in section 5 and eq. (3.2) of 
[52] provides a non-trivial check of the results in [50, 58]. 

Before concluding we want to comment on interesting questions regarding the 
construction and the universality of the e-finite basis. 

The construction of the e-finite basis as discussed in section 3 requires that 
the reduction of the given set of Feynman integrals to some basis of master in- 
tegrals has already been performed. For purely technical reasons this reduction 
procedure is essentially more complicated and time consuming than finding the 
e-finite basis. A significant part of the complexity arises from the possible ap- 
pearance of spurious poles since, as a consequence, the exact dependence on d of 
all coefficient functions has to be kept through the whole reduction procedure. 
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In a modified approach one could try to implement the idea of the e-finite 
basis directly in the Laporta Algorithm, by e.g. properly tuning specific weights 
in the algorithm. This would lead to a new criterion instead of the commonly 
used one, which is based, roughly speaking, on the number of lines, dots and the 
power of irreducible numerators. This new criterion could, in addition, depend 
on the appearance of {d — 4) in the coefficients of the integrals, something like 
"never solve an equation with respect to an integral, which has a factor of {d — 4) 
in its coefficient" . As a result no spurious pole would appear in any step, which 
would allow to expand the coefficient functions in all equations around the point 
d = 4 up to some fixed power in e — {A — d)/2 and would, in general, lead to a 
significant gain in performance of the Laporta Algorithm. 

On the other hand, it might happen that the IBP identities "globally" propor- 
tional to some positive power of {d — 4) play an important role and that without 
their use no complete reduction is possible. The problem is open and certainly 
deserves further investigation. 

The second interesting problem is whether the e-finiteness of a basis is pre- 
served if one considers a different physical problem whose solution is expressible 
in terms of the "old" master integrals. On general grounds we can not exclude 
that for an "old" master, say, Tj the maximal power of the spurious pole would 
e.g. exceed the one displayed on figure 2. In such the basis given in section 

4.1 docs not have coefficient functions finite in the limit e ^ and we would 
have to construct a new, more general e-finite basis following the prescription of 
section 3. However, we do think that such a situation is unlikely and that there- 
fore the in section 4.1 constructed e-finite basis is a universal one for all four-loop 
problems reducible to the evaluation of QED hke massive tadpoles. The following 
observations support our point of view: 

• As stated in section 4 we have only used the reduction and integrals ap- 
pearing in the low energy expansion of the four-loop vacuum polarization 
function (see [34]) for constructing the e-finite basis. In fact, the evaluation 
of decoupling function Q for the QCD quark gluon coupling constant in 
four-loop approximation requires the reduction of additional more compli- 
cated integrals to masters. We have checked that the basis given in figure 3 
remains e-finite for all these extra integrals. 

• If our hypothesis on the universality property of the constructed e-finite 
basis is wrong it also implies that using the method described in section 
4.3 one should be able to find more terms in the e-expansions of the "old" 
and of the e-finite master integrals analytically. This extension of the equa- 
tions (24-32,34-45,51) without using more sophisticated techniques seems 
not very probable to us. 
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Conclusions 

As a conclusion we want to point out that the very idea of the e-finite basis, 
namely the procedure of its systematic construction as well as the application to 
finding direct analytic information of e-finite contributions to non-e-finite master 
integrals are all new. We beheve that choosing the e-finite basis should be es- 
pecially useful and convenient in all cases when numerical evaluation of master 
integrals is necessary. 

Our results are also available in computer readable form under the URL 
http : / /www-ttp . physik . uni-karlsruhe . de/Progdata/ttp06/ ttp06-02/. 
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Old masters in terms of e-finite basis 



52 — 

2 (rf - 2)2 (1298 - 953 + 233 ^2 - 19 (^3) 

^ {d-3){3d- 10) (3 - 8) (73 - 34 d + 4 d^) 
4(d-5) (d-4) (4d-17) 
(2d-7) (73-34d + 4d2) 

16 (2 d - 7) (2 d - 5) (209 - 96 d + 11 d^) 
^ 3(d-2) (3rf- 11) (73 - 34ci + 4d2) 

32 (d - 4) (rf - 3)2 (2 - 7) (149 - 69d + 8£) 
^ (d - 2) (3 ci - 11) (3 - 10) (3 - 8) (73 - 34 ci + 4 ci^) ' 



32 (ci - 5) (ci - 4)2 (rf - 3) 
(d-2)(3d-10)(3rf- 8)(73 - 34(i + 4(i2) ^2 

256 (ci-5) (d-4)^ f 

3{d-2) {3d- 10) (3(i- 8)(73- 34d + 4d2) 

16 (d- 4)3 (767 - 351^ + 40^2) 
3(d-2) (3(^-10) (3rf- 8) (73- 34(i + 4d2) ei 
128(r/-5)2(r/- 1)^^ . 



-^62 



(d - 2) (3 d - 10) (3 d - 8) (73 - 34 d + 4 ^2) 

32 (fi- 4)^ (109-51^ + 6^2) f 

{d-2){3d- 10) (3 d - 8) (73 - 34 d + 4 d^) 
48(rf-4)3(149- 69(i + 8ci2) 

(d - 2) (3 ci - 10) (3 - 8) (73 - 34 d + 4 d"^) 
16(ci-4)3(673 - 297rf + 32rf2) 

3 (d - 2) (3 - 10) (3 d - 8) (73 - 34 d + 4 rf2) ' 

128 {d-5) {d-Af{3d-13) 
3(d-2)(3d- 10) (3 - 8) (73 - 34 d + 4 d'^) 



'-91 
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{d - 2f (21458 - 20723 d + 7488 (f - 1200 + 72 d^) 
2(c^-3) (2ci-7)(2d-5) (3ci-10) (3d -8) (73- 34^ + 4^2) 
8 (d - 4)2 (rf - 2) (25 - 22rf + 4rf2) 
^ (2(i-7)2(2d- 5) (73- 34ci+4ci2) 
32(d-4) (lOrf-33) 



+ 
+ 



3(3(i-ll)(73- 34rf + 4rf2) 

256 (rf-4)2(rf-3)3 



63 



{2d -5) {3d- 11) (3 - 10) (3 - 8) (73 - 34 d + 4 d'^) 

256(d-5)(d-4)3(d-3) ^ 

{2d -7) {2d -5) {3d- 10) (3 (/ - 8) (73 - 34 + 4 ^2) ^2 

512 (r/- 1)' . 

3(2(i-5)(3ci-10)(3ci- 8)(73- 34(i + 4(i2) 

128(^-4)^(11^-37) J 

3(2ci-7)(2ci-5)(3d-10)(3d- 8)(73- 34d + 4(i2) ei 

+ 



+ 



256 {d -5){d- 4)3 



(2(i-5) (3(i-10) (3rf- 8)(73 - 34(i + 4(i2) 



128(^-4)^ f 

----- ----- - - . -^64 



{2d -7) {2d -5) {3d- 10) (3 - 8) (73 - 34 d + 4 d'^) 

384 {d - 4)^ {d - 3) 

{2d - 7) {2d - 5) {3d - 10) {3d - 8) {73 - 3Ad + Ad'^) 

128 (d - 4)^ (25 d - 89) . 



3 (2 d - 7) (2 d - 5) (3 d - 10) (3 d - 8) (73 - 34 d + 4 d"^) 

256 (d- 4)=^ (3d -13) 

3(2d-5)(3d-10)(3d- 8)(73- 34d + 4d2) 
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^64 — 

3(d-4) {d-2) {3d' 8) (82-48d + 7d^) 

64 (rf-3)3(2rf-7)2 ^ 
(2d -5) (3d- 10)2 (3d -8) 

4 (d - 4) (d - 3)3 
3(3d-10) ^ 

4(d-3) 
9 (d - 4)2 (3 d - 11) (3 d - 10) . 
8(d-3)3(2d-7) 



20 



3(d-2)^(11278- 11551 d + 4411 d^- 745 + 47 d^) 

^ 16 {d - 3)3 (2 - 7) (3 - 10) (73 - 34 + 4 d^) 
3 (d - 5) (c^ - 4) {d-2){3d- 8) (4 d - 17) (5 ci - 16) 

^ 8 (rf - 3)2 (2 - 7)2 (73 - 34 d + 4 rf2) 

(2 d - 5) (3 - 8) (21406 - 22007 d + 8444 d^ - 1434 rf^ + 91 d'^) 

2 (rf - 4) (rf - 3)2 (3 - 1 1) (73 - 34 + 4 rf2) 
3 (17566 - 18119 d + 6980 d'^ - 1191 d^ + 76 d^) 



+ 
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(3d - 11) {3d- 10) (73-34d + 4c^2) 

3(d-5)(d-4)2(5d-16) f 

{d- 3) {2d -7) {3d - 10) (73 - 34 d + 4 d'^) 

((i-3)2(2d-7)(3(i- 10)(73-34d + 4d2) 54 

{d - 4)2 (73178 - 75093 d + 28748 d^ - 4869 ci^ + 308 d^) . 



2((i-3)2(2rf- 7) (3d - 10) (73 - 34d + 4(i2) 
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12 (ci-5)2(ci-4)3(5ci- 16) 



{d - 3)2 {2d -7) {3d- 10) (73 - 34 d + 4 ^2) ^62 
3{d- 4)2 (15006 - 15527 d + 6004 ^2 _ 1029 d^ + 66 d^) 

(d- 3)2 (2d -7) (3d - 10) (73- 34d + 4d2) ^ 
9 (d - 4)2 (17566 - 18119 d + 6980 d2 - 1191 d^ + 76 d^) f 

2 (d - 3)2 (2 d - 7) (3 d - 10) (73 - 34 d + 4 d^) 
(d - 4)2 (18982 - 18219 d + 6388 d2 - 963 d^ + 52 d^) . 
2 (d- 3)2 (2d -7) (3d - 10) (73- 34d + 4d2) 
4(d-5) (d-4)3(3d- 13)(5d- 16) 
(d - 3)2 (2 d - 7) (3 d - 10) (73 - 34 d + 4 d2) ' 
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4(ci-3)3(2c^- 7)(73- 34d + 4d2) 

(d - 2) (3 - 8) (499 - 409 + 111 rf^ _ lo ^3) 

4(d-3)2(2rf- 7) (73 - 34rf + 4rf2) 
2 (2 d - 5) (3 - 10) (3 - 8) (63 - 28 + 3 d^) 

3((i-3)2(3d- 11) (73- 34 (i + 4(^2) 
4(c^-4)2(4rf- 19) 



63 



(3d - 11) (73 - 34 ci + 4 (^2) 

4(d-5)(d-4)2 . 
(d-3)(2d- 7)(73- 34d + 4d2) 52 

16(d-4)3(33- 16d + 2d2) 
3 (d- 3)2 (2d -7) (73-34d + 4d2) 54 
2 (d -4)3 (255- 125 (i+ 16 (i^) . 
3 (d- 3)2 (2d - 7) (73- 34d + 4(i2) ei 

4(d-5)(d-4)2 . 
(d-3)2(73 -34d + 4d2) 62 

4(rf_4)4(2rf-9) 
(d-3)2(2d- 7)(73 -34d + 4d2) 64 
6 (d- 4)^4^-19) ; 
(d-3)2(2d- 7)(73-34d + 4d2) 

2(d-4)3(276- 145d + 20d2) . 
3 (d- 3)2 (2d -7) (73-34d + 4d2) 
8 (d- 4)^ (3d - 13) (33- 16d + 2d^) 
3 (d- 3)2 (2d -7) (73-34d + 4d2) ' 
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63 



{d - 2f (1789 - 2111 d + 912 - 172 + 12 d^) 

8 (ci - 3)2 (2 - 7) (3 - 10)2 (73 - 34 d + 4 d^) 
{d - 4) (d - 2) (3 ci - 8) (914 - 525 + 90 rf2 _ 4 d^) 

16 (d - 3) (2 - 7)2 (3 ci - 10) (73 - 34 ci + 4 d^) 
(2 d - 5) (3 - 8) (9772 - 7470 d + 1165 d'^ + 399 rf^' - 142 rf^ + 12 d^) 

6(rf-4)2(rf-3) (3 rf- 11) (3d - 10) (73 - 34(i + 4(i2) 
(d - 3)2 (12628 - 13328 d + 5255 d"^ - 918 d^ + 60 c^^) 
{d-A) (3d -11) (3rf- 10)2 (73- 34d + 4rf2) 
S{d-h){d-Af f 
(2d-7)(3d- 10)2(73- 34d + 4d2) 52 

16(^-4)^ / 

3(d-3)(3d- 10)2(73- 34d + 4d2) 54 

4 (d- 4)^ (lid -37) f 

3{d-3) (2rf-7)(3rf- 10)2 (73 -34d + 4ci2) ei 

8{d-5){d-4f f 
(rf-3)(3rf- 10)2(73 - 34rf + 4rf2) 62 

2 (d - 4) (35898 - 47927 d + 25554 d'^ - 6802 rf^ + 904 rf^ - 48 d^) . 

(rf-3) (2 a! -7) (3rf- 10)2 (73 - 34 ci + 4 rf2) ' 

3 {d-A) (12628 - 13328 d + 5255 (^2 - 918 d^ + 60 d^) f 

2 (2d -7) (3d - 10)2(73- 34d + 4d2) 

4(d-4)^(25d-89) f 

3 (d - 3) (2 d - 7) (3 d - 10)2 (73 - 34 d + 4 d'^) 

8(d-4)^(3d-13) 
3(d-3)(3d- 10)2(73- 34d + 4d2) ' 



64 
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Tn = (67) 
{d - 2)3 (5925 - 6293 d + 2499 d^ - 440 d^ + 29 d^) 

4{d- 3)2 (2 - 7) (3 - 10)2 (73 - 34 + 4 d^) 
{d-2){3d- 8) (9053 - 9587 d + 3799 d"^ - 668 d^ + 44 d^) 

4 (rf - 3) (2 c/ - 7)2 (3 - 10) (73 - 34 + 4 d^) 
2 (2 - 5) (3 - 8) (26546 - 35605 d + 19114 d^ - 5132 ci^ + 689 rf^ - 37 rf^) 

3(rf-4)2(d-3) (3ci- 11) (3rf- 10) (73 - 34 d + 4 ^2) ' 
4 (rf - 3) (19090 - 25909 d + 14064 d^ - 3815 d^ + 517 rf^ - 28 d^) 
(c/ - 4) (3 - 11) (3 c/- 10)2 (73 - 34 + 4 d^) 
4 (rf - 5) (c/- 4)2(5^-19) J 
(2 - 7) (3 - 10)2 (73 - 34 + 4 ^2) 52 
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■T, 
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16 (rf- 4)3 (554 -435 c/ +114^2 -10^3) 
^ 3 (c/ - 3) (2 ci - 7) (3 ci - 10)2 (73 - 34 c/ + 4 c/2) ^4 
2[d- 4)3 (7327 - 5682 d + 1473 c/2 - 128 d^) . 
^ 3{d-3){2d-7){3d- 10)2 (73 - 34 c/ + 4 c/2) 
4(c/-5) (ci-4)2 (54 - 29c/ + 4c/2) . 
(c/ - 3) (3 c/ - 10)2 (73 - 34 c/ + 4 c/2) ^2 

4 (c/ - 4) (21906 - 29301 d + 15680 c/2 - 4195 + 561 c/^ - 30 d^) f 
^ {d-3){2d~l){3d- 10)2 (73 - 34 c/ + 4 c/2) ^'^^ 

6 (c/ - 4) (19090 - 25909 c/ + 14064 d'^ - 3815 d^ + 517 c/^ - 28 d^) , 
^ (c/ - 3) (2 c/ - 7) (3 c/ - 10)2 (73 - 34 ci + 4 d^) ' 

2 (c/ - 4) (50762 - 63841 d + 32144 c/2 - 8107 c/^ + 1025 c/^ - 52 d^) . 
3{d-3){2d-l){3d- 10)2 (73 - 34 c/ + 4 d'^) 

8 (ci - 4)2 (3 c/- 13) (554 - 435 c/ + 114 c/^ - 10 d^) 
^ 3{d-3){2d-7){3d- 10)2 (73 - 34 c/ + 4 c/2) ' 

B Proceeding in case of a "divergent" cut 

In this section we show an alternative approach of calculating diagram Tg^. If we 
perform the cut of the diagram in the following way 

S^-(P)= P (68) 



one obtains a divergent self energy function S-^52(p). This results from the fact 
that the diagram in equation (68) consists of 



X ^^ (69) 

^ v ' 

finite =i+finite 
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a finite two-loop multiplied by a divergent one-loop diagram. If one identifies 
such a sub- divergency one is able to construct the combination 



X 



- +- P 



(70) 



eS^52 (p)=fimte 



of a finite self energy function S^sa (p) and a diagram of lower loop order which 
can be calculated analytically. In general one might have to treat several such 
sub- divergencies in the similar way. In a next step E^52 (p) 



E^^(p) 



1 

- P 
e 



(71) 



is evaluated to obtain T^2[n/m] according to the proceeding described in section 
4.2. Thus we are again able to extract the pole structure of T^2[n/m] analytically 
and the constant part of the e-expansion numerically. Finally one obtains the 



expansion for the diagram Tc 



52 



52 [n/m] 



/ 



52 [n/m] 



+ - 




(72) 



by adding the analytical result of the three-loop tadpole diagram. We find full 
agreement between this approach and the results discussed in section 4.2. 
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